Suppression of non-adiabatic losses of molecules from chip-based microtraps 
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Polar molecules in selected quantum states can be guided, decelerated, and trapped using electric 
fields created by microstructured electrodes on a chip. Here we explore how non-adiabatic transitions 
between levels in which the molecules are trapped and levels in which the molecules are not trapped 
can be suppressed. We use ^^CO and ^^CO [a ^111 ,v = 0) molecules, prepared in the upper A-doublet 
component of the J = 1 rotational level, and study the trap loss as a function of an offset magnetic 
field. The experimentally observed suppression (enhancement) of the non-adiabatic transitions for 
^^CO ('^^CO) with increasing magnetic field is quantitatively explained. 

PACS numbers: 37.10.Pq, 31.50.Gh, 37.10.Mn 



(N 

Oh. 



: 

o ■ 

^ : 

Oh. 



> 

(N 
O 



I. INTRODUCTION 



The manipulation and control of polar molecules above 
a chip using electric fields produced by microstructured 
electrodes on the chip surface is a fascinating new re- 
search field Miniaturization of the electric field struc- 
tures enables the creation of large field gradients, i.e., 
large forces and tight potential wells for polar molecules. 
A fundamental assumption that is made when consider- 
ing the force imposed on the molecules is that their po- 
tential energy only depends on the electric field strength. 
This is usually a good assumption, since the molecules 
will reorient themselves and follow the new quantization 
axis when the field changes direction and their potential 
energy will change smoothly when the strength of the 
field changes. This approximation can break down, how- 
ever, when the quantum state that is used for manipula- 
tion couples to another quantum state that is very close 
in energy. If the energy of the quantum state changes 
at a rate that is fast compared to the energetic split- 
ting, transitions between these states are likely to oc- 
cur. For trapped molecules in so-called low-field-seeking 
states in a static electric potential, such transitions are 
particularly disastrous when they end up in high-field- 
seeking states or in states that are only weakly infiu- 
enced by the electric fields, as this results in a loss of 
the molecules from the trap. This effect has been in- 
vestigated previously for ammonia molecules in a loffe- 
Pritchard type electrostatic trap with a variable field 
minimum. In this macroscopic electrostatic trap, losses 
due to non-adiabatic transitions were observed on a sec- 
ond time scale when the electric field at the center of the 
trap was zero; with a non-zero electric field minimum 
at the center of the trap, these losses could be avoided 
Non-adiabatic transitions have recently also been 
investigated in a "conventional", i.e. macroscopic. Stark 
decelerator in which electric fields are rapidly switched 
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between two different configurations. There, these tran- 
sitions have been found to lead to significant losses of 
molecules when they are in low electric fields Q . Similar 
trap losses will be much more pronounced on a microchip, 
where the length scales are much shorter and where the 
electric field vectors change much faster. For atoms in 
a three-dimensional magnetic quadrupole trap, the trap 
losses due to spin flip (or Majorana) transitions has been 
shown to be inversely proportional to the square of the 
diameter of the atom cloud On atom-chips, where 
paramagnetic atoms are manipulated above a surface us- 
ing magnetic flelds produced by current carrying wires, 
trap losses due to Majorana transitions are therefore well- 
known but can be conveniently prevented by using an 
offset magnetic field 01 ■ Due to the geometry of the 
molecule chip, however, applying a static offset electric 
field is not possible, and other solutions must be sought. 

We have recently demonstrated that metastable CO 
molecules, laser-prepared in the upper A-doublet com- 
ponent of the J = 1 level of the a'^IIi,'!; = state can 
be guided, decelerated and trapped on a chip. In these 
experiments, non-adiabatic losses have been observed for 
12(-j16q Jj-^ ^jjjg j^-^ogi; abundant carbon monoxide isotopo- 
logue, the level that is low-field-seeking becomes degen- 
erate with a level that is only weakly influenced by an 
electric fleld when the electric held strength goes to zero. 
Every time that the trapped molecules pass near the zero 
fleld region at the center of a micro-trap, they can make 
a transition between these levels and thereby be lost from 
the trap. This degeneracy is lifted in ^-^C^^O due to the 
hyperfine splitting (the ^^C nucleus has a nuclear spin 
|/| = 1/2), and the low-field-seeking levels never come 
closer than 50 MHz to the non-trappable levels. There- 
fore, changing from ^^C^^O to ^^C^^O (referred to as 
^^CO and ^'^CO from now on) in the experiment greatly 
improves the efficiency with which the molecules can be 
guided and decelerated over the surface and enables trap- 
ping of the latter molecules in stationary traps on the 

chip [3. 

Although it is evident that the 50 MHz splitting be- 
tween the low-field-seeking and non-trapped levels in 
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^^CO is beneficial, it is not a priori clear whether a 
smaller splitting would already be sufficient or if a still 
larger splitting would actually be needed to prevent all 
losses. While the hyperfine splitting in ^'^CO cannot be 
varied, the degeneracy can be lifted by a variable amount 
in the normal ^^CO isotopologue by using a magnetic 
field. If a magnetic field is applied in addition to the 
electric field, a splitting can be induced between the low- 
field-seeking and non-trappable levels of ^^CO that de- 
pends on the strength of the applied magnetic field; in 
^■^CO, a magnetic field will actually decrease the splitting 
between the low-field-seeking and non-trappable levels. 

In this paper, we present measurements of the effi- 
ciency with which CO molecules are transported over the 
chip — while they are confined in electric field minima 
that are traveling at a constant velocity — as a function 
of magnetic field strength. It is observed that, in the 
case of ^^CO, the losses due to non-adiabatic transitions 
can be completely suppressed in sufficiently high mag- 
netic fields; for ^•^CO on the other hand, the trap losses 
increase with increasing magnetic field. A theoretical 
model that can quantitatively explain these observations 
is also presented. Although the chip is ultimately used 
to decelerate molecules to a standstill, measuring the ef- 
ficiency with which molecules are guided at a constant 
velocity provides a detailed insight into the underlying 
trap-loss mechanism. Limiting the experiments to con- 
stant velocity guiding also makes them more tractable: 
bringing molecules to a standstill and subsequently de- 
tecting them has thus far required five separate phases of 
acceleration, which greatly complicates efforts to under- 
stand the details of the loss mechanism using numerical 
calculations While deceleration at a constant rate 
to a non-zero final velocity is possible, the measurable 
signal in such experiments is significantly lower than in 
constant velocity guiding. Despite the fact that we only 
measure at constant velocity, we nonetheless apply the 
model to examine the non-adiabatic losses expected dur- 
ing deceleration. This work thereby furthers the goal of 
extending trapping on the molecule chip to a wider range 
of molecules. 



II. EXPERIMENTAL SETUP 

A scheme of the experimental setup is shown in Fig- 
ure [T] A mixture of 20% CO in krypton is expanded into 
vacuum from a pulsed valve (General Valve, Series 99), 
cooled to a temperature of 140 K. In this way, a molecular 
beam with a mean velocity of 300 m/s and with a full- 
width-half-maximum spread of the velocity distribution 
of approximately 50 m/s is produced. This beam passes 
through two 1 mm diameter skimmers and two differen- 
tial pumping stages (the valve and first skimmer are not 
shown in the figure) before entering the chamber in which 
the molecule chip is mounted. Just in front of the second 
skimmer, the ground state CO molecules are excited to 
the upper A-doublet component of the J = 1 level of the 
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FIG. 1: Schematic of the experimental setup. The laser pre- 
pared metastable CO molecules are guided over the chip in 
tubular electric field minima that move at a constant speed of 
300 m/s. Above the chip, the molecules are exposed to a time- 
dependent electric field as well as to a constant magnetic field 
whose direction is always perpendicular to the electric field. 
The guided CO molecules are ejected from the chip, ionized 
via laser induced resonance enhanced multi photon ionization 
(REMPI), and the CO"*" ions are mass-selectively detected on 
a micro-channel plate (MCP) detector. The coordinate sys- 
tem used in the further description is explicitly indicated. 



metastable a'^IIi, v = state, using narrow-band pulsed 
laser radiation around 206 nm (1 mJ in a 5 ns pulse 
with a bandwidth of about 150 MHz). The metastable 
CO molecules are subsequently guided in traveling po- 
tential wells that move at a constant speed of 300 m/s 
parallel to the surface of the molecule chip. A uniform 
magnetic field is applied to the region around the chip 
using a pair of 30 cm diameter planar coils separated 
by 23 cm (not shown in the figure). The coils are ori- 
ented such that the magnetic field is parallel to the long 
axis of the chip electrodes, i.e. along the z-axis, ensuring 
that the magnetic field is always perpendicular to the 
electric field {vide infra). The CO molecules that have 
been stably transported over the chip will pass through 
the 50 Jim high exit slit and enter the ionization detec- 
tion region a short distance further downstream. There, 
the metastable CO molecules are resonantly excited to 



selected rotational levels in the 6 



state us- 



ing pulsed laser radiation at 283 nm (4 mJ in a 5 ns 
pulse with a 0.2 cm~^ bandwidth). A second photon from 
the same laser ionizes the molecules and the parent ions 
are mass-selectively detected in a compact linear time-of- 
flight setup using a micro-channel plate (MCP) detector. 
This detection scheme has been implemented in addition 
to the Auger detection scheme that we have used in ear- 
lier studies P, 0, 01 as it is more versatile and can also be 
applied to detect other molecules. In addition, the de- 
tection sensitivity of the ion detector is less affected by 
the magnetic field than that of the Auger detector. 

The molecule chip and its operation principle have 
been described in detail before 0| , and only the fea- 
tures that are essential for understanding of the present 
experiment are discussed here. The active area of the 
chip consists of an array of 1254 equidistant electrodes. 
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FIG. 2: Edge-on view (in the +2-direction) of the molecule 
chip with calculated contour lines of equal electric field 
strength above the chip, displayed at intervals of 0.45 kV/cm. 
The three panels correspond to three different times in the 
harmonic waveform cycle. The position of the electrodes is 
indicated at the bottom of each panel and the instantaneous 
values of the applied potentials (in volts) are given. The time 
difference between adjacent panels is j^, where u is the fre- 
quency of the harmonic waveforms. From the inset in each 
panel, it is seen that the electric field at each position relative 
to the center of the trap rotates clockwise over an angle of j 
between adjacent panels. 



each 10 pm wide and 4 mm long, with a center-to-center 
distance of 40 pm. An edge-on view of the chip electrodes 
(with the 4 mm dimension of the electrodes perpendic- 
ular to plane of the figure) is shown in Figure [21 The 
potential (in volts) applied to an electrode at a given 
moment in time is indicated directly above the electrode; 
these six potentials are repeated periodically on the elec- 
trodes on either side of those drawn here. Because the 
electrodes are much longer than the period length of the 
array, the electric field distribution can be regarded as 
two-dimensional, i.e. the component of the electric field 
along the z-axis can be neglected. This is of importance 
for the present experiments, because only in this case the 
applied magnetic field is always perpendicular to the elec- 
tric field. The calculated contour lines of equal electric 
field strength in the free space above the chip show elec- 
tric field minima that are separated by 120 pm, i.e., there 
are two electric field minima per period, centered about 
25 pm above the surface of the chip. By applying sinu- 
soidal waveforms with a frequency v to the electrodes, 
these minima can be translated parallel to the surface 
with a speed given hy v ~ 120 pm • v. When these wave- 
forms are perfectly harmonic and have the correct ampli- 



tude, offset, and phase, the minima move with a constant 
velocity at a constant height above the surface, and the 
shape of the field strength distribution does not change 
in time. Because an electric field strength minimum acts 
as a trap for molecules in low-field-seeking states, these 
fields act as tubular moving traps that can be used to 
guide the molecules over the surface of the chip. The 
tubular traps are closed at the end by the fringe fields 
caused by the neighboring electrodes. Near the ends of 
the about 4 mm long traps, the electric field will neces- 
sarily have a component along the z-axis, i.e. parallel to 
the applied magnetic field. In the present study, where 
the molecules are guided at 300 m/s over the chip and 
are therefore on the chip for less than 200 ps, these end- 
effects are neglected. 

The region near an electric field minimum at {x,y) = 
{xo,yo) is a quadrupole. with an electric potential given 
by 

y = |r2cos(2^-</>o) (1) 

where x — xo = r cos (p and y — yo = r sin (p. In the current 
experiment, sinusoidal waveforms with an amplitude of 
180 V are applied to the electrodes, yielding a value of 
a = 0.054 V/pm^ 0|. The resulting electric field is given 
by 

E = — ar(cos(0o ~ 0)^ + sin((/io — </')y)- (2) 

The strength of the electric field, \E\ = ar, depends only 
on the r coordinate, but the direction of the field vector, 
(f)' — (j)o — (j) + TT , depends on the coordinate and on the 
phase factor c/jq . While the direction of the field vector 
changes as a result of the motion of the molecule in the 
quadrupole field (and thus changing (f>), the direction of 
the field at any given position relative to the minimum 
also rotates when the minimum is translated over the 
chip. It is seen from the electric field vectors shown in 
the insets of Figure [2] that the frequency of this rotation 
is 1.5 times the frequency of the applied waveforms and 
that the direction of the rotation is clockwise, i.e. the 
rotation vector points along the positive z-axis and (po 
increases linearly in time. To guide the molecules over the 
chip at 300 m/s, harmonic waveforms with a frequency 
J/ of 2.5 MHz must be applied, resulting in a rotation 
frequency of 3.75 MHz. 

For CO molecules in the low-field-seeking component 
of the J = 1 level, the depth of the tubular traps above 
the chip is about 60 mK. This implies that CO molecules 
with a speed of up to 6 m/s relative to the center of the 
trap can be captured. The oscillation frequency of the 
molecules in the radial direction of the tubular traps is 
in the 100-250 kHz range. These parameters are of im- 
portance for the non-adiabatic transitions as these deter- 
mine with which velocity, and how often per second the 
CO molecules pass by the zero-field region of the traps. 

It turns out that, in the actual experiment, the tubu- 
lar traps do not move perfectly smoothly over the chip. 
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FIG. 3: Calculated position of the minimum of an electric field 
trap with respect to an ideal trap moving at constant veloc- 
ity during an 800 ns time-interval together with the resulting 
velocity in the xy plane, using the measured "standard wave- 
forms" (upper row) and "improved waveforms" (lower row) as 
input. 



Due to imperfections in the amplitude, offset and phase 
of the waveforms that we have used, the tabular traps 
are jittering at rather high velocities. The motion of the 
center of the traps relative to the ideal, constant velocity 
motion can be calculated by measuring the real wave- 
forms applied to the chip and using these to compute 
the position of the minimum for each point in time. The 
upper row of Figure [3] shows the motion of a minimum 
when the waveforms that we will refer to as the "standard 
waveforms" have been used. The range of this motion ex- 
tends over ±2 p.m in the x and y coordinates. Although 
this is considerably smaller than the size of the trap- 
ping region, this motion significantly enlarges the effec- 
tive region in which non-adiabatic transitions can occur. 
Moreover, the entire path is traced out periodically every 
T ^ 2/iy ~ 800 ns; because the motion occurs on such 
a short time scale, the speed with which the trap center 
moves, and therefore the relative speed with which the 
molecules encounter the trap center, can be as high as 
100 m/s. To improve the waveforms, we inserted an LC 
filter in the output stage of the amplifiers, thereby reduc- 
ing the harmonic distortion. The resulting motion using 
these "improved waveforms" is shown in the lower row of 
Figure [31 It is seen that not only the range of the motion 
is now contracted but that also the speed with which the 
trap center jitters is reduced by about a factor of two. 



III. THEORETICAL MODEL 

A. Eigenenergies in combined fields 

In order to describe the non-adiabatic transitions in 
CO, we must first derive the energy levels of CO in com- 
bined electric and magnetic fields. For this, the field-free 
Hamiltonian is expanded with Stark and Zeeman contri- 
butions, i.e., 

H = i?A,hfs + Hs + Hz- (3) 

Here, Hh,\\fs describes the A-doubling of the a '^IIi, w = 0, 
,7 = 1 level for either ^^CO or ^"^CO and also includes the 
hyperfine splitting of each A component into F = 1/2 and 
F = 3/2 hyperfine sublevels for "CO; Hs = - 'Ue ■ E is 
the Stark interaction Hamiltonian and Hz — — M ' B is 
the Zeeman interaction Hamiltonian, where JIe and /i are 
the electric and magnetic dipole moment operators, re- 
spectively, and E and B are the (time-dependent) electric 
and magnetic field vectors. 

The spectroscopic parameters of the a'^Hi, u = state 
of CO that are used in the field-free Hamiltonian are 
given elsewhere (sl-fisj. The A-doublet splitting between 
the positive parity component (upper) and the negative 
parity component (lower) of the J = 1 level is about 
400 MHz while the hyperfine splitting of each parity level 
of "CO into F = 1/2 (lower) and F = 3/2 (upper) 
sublevels is about one order of magnitude smaller. The 
body-fixed electric dipole moment ^le = ImbI in the elec- 
tronically excited metastable state is 1.3745 Debye for 
both I^CO and "CO The magnetic moment of 

the molecule can be expressed as ^ = —/is [gL-L+gs-S), 

where /xb is the Bohr magneton, L is the electron orbital 

angular momentum operator, S is the electron spin op- 
erator, and where the magnetic g-factors are fixed at the 
values of the bare electron, g^ = \.Q and gs = 2.0023. 

A detailed description of the formalism used to calcu- 
late the eigenenergies of the various components of the 
J = 1 level in the a^Hi, v — Q state of both ^^CO and 
^■^CO in combined electric and magnetic fields is pre- 
sented in the Appendix. The formalism has been set 
up for mutually orthogonal static electric and magnetic 
fields. As will be discussed below, this is also adequate 
to treat the actual situation, in which the electric field 
rotates with a constant frequency in a plane perpendic- 
ular to the magnetic field. In Figure 0] we only show 
the outcome of these calculations in the form of plots of 
the energy levels for the upper A-doublet component of 
^^CO and ^'^CO as a function of electric field strength 
in the absence (left column) or in the presence of a 
50 Gauss magnetic field (right column). At low electric 
field strengths, the Stark shift is quadratic, but as the 
electric field strength increases and \Hs\ becomes much 
larger than |i?A.hfs| and \Hz\, the Stark energy shows 
a linear dependence on the electric field strength, i.e. 
Ai?s oc —ilMj fj.EE, and the product of the projection of 
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FIG. 4: Energy of the upper (positive parity) A-doublet com- 
ponent of the J = 1 level of ^^CO (upper panels) and ^^CO 
(lower panels) in the electronically excited a ^111 , « = state 
as a function of the electric field strength, without (left col- 
umn) and with (right column) an offset magnetic field. The 
strength of the magnetic field is 50 Gauss and its direction 
is perpendicular to the direction of the electric field. The 
labeling of the energy levels is explained in the text. 



the electronic angular momentum along the internuclear 
axis 17 and the projection of the angular momentum 

J on the electric field vector E becomes an approximately 
good quantum number. If a weak magnetic field such as 
that present in the experiment is applied in the absence 
of an electric field, each of the zero-field levels splits into 
(2 J-l-1) (or in (2_F-|-1)) separate levels based on their 
(or MJ?) quantum number. The Zeeman energy is linear 
in magnetic field and can be calculated using first order 
perturbation theory as l^Ez ~ ficsMB where M = Mf 
(or Mp) and fi^s ^ f^B describes the effective magnetic 
moment of a particular parity and, in the case of ^^CO, 
particular F component. When a strong electric field 
is applied to the molecule in addition to the magnetic 
field, such that \Hs\ ^ |7?A^hfs|, the energy levels 

can again be characterized with the approximately good 
quantum number U,Mf and show a linear Stark shift. 

The behavior of the eigenenergies in static electric and 
magnetic fields is quite instructive in describing the non- 
adiabatic transitions — and thus losses — of molecules 
from the low-field-seeking states {^Mf = —1) to non- 



trappable {Mj = 0) states. Near the edge of the trap, 
where the electric field is as high as 4.2 kV/cm, the 
Stark effect provides an energy gap of about J/trap = 
60 mK ~ 1.3 GHz between low-field-seeking states and 
non-trappable states. Since this is much larger than 
both the frequency of the motion of the molecules in the 
traps and the frequency of the applied waveforms, non- 
adiabatic losses will not occur near the edge of the trap. 
In the vicinity of the trap center, however, this argu- 
ment no longer holds, and the eigenfunctions can change 
from A/|'-type wavefunctions to Mf -type (or A/J?-type) 
at a rate faster than the energy gap iJgap in that region. 
In the case of a ^^CO molecule, the energy gap at the 
center of the trap goes to zero in the absence of a mag- 
netic field. If a molecule in a low-field-seeking state flies 
near the trap center with a velocity v high enough, or 
with a distance of closest approach b small enough, that 
the corresponding interaction time with the trap center 
T = b/v no longer fulfills the adiabaticity condition, i.e. 
when the condition iJgap ^ h/r no longer holds, then 
the probability of transitions to non-trappable states can 
become significant. In the absence of a magnetic field. 
^''CO molecules are much safer from such non-adiabatic 
losses due to the energy gap of 50 MHz between the 
F ^ 3/2 level (which becomes low- field- seeking in an 
electric field) and the F = 1/2 level (which correlates 
with non-trappable Mf = states). 



B. Calculating rates for non-adiabatic transitions 

The basic idea underlying the calculation of the non- 
adiabatic losses for the molecules in low-field-seeking 
states is that, since the transitions to non-trappable 
states happen primarily as the molecules pass the zero 
field region at the center of a microtrap, the overall loss 
probability can be estimated by first calculating the loss 
probability in a single pass. For simplicity, the trajec- 
tory of the CO molecules is assumed to have a con- 
stant velocity: this is reasonable, since the forces on the 
molecules approach zero at the center of the trap due to 
the A-doubling. The transition probability P,;.j (v, B, b) of 
a molecule making a transition from a low-field-seeking 
state i to a non-trappable state j for a single pass by the 
trap center depends on the speed of the molecule relative 
to the center of the trap v, on the strength of the mag- 
netic field B, and on the distance of closest approach b. 
Due to the rotation of the electric field, there is a differ- 
ence between positive and negative values of b, which is 
why we do not refer to it as an "impact parameter" here. 

To calculate the probability Pi,j{v,B,b), the time- 
dependent Schrodinger equation, which describes the 
evolution of the quantum states, must be solved: 



ih— = H^. 
at 



(4) 



This equation is solved numerically as an initial value 
problem on a set of coupled first-order ordinary differ- 
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ential equations using the basis vectors given in the Ap- 
pendix. The Hamihonian H depends on the electric field 
vector E, which, for a molecule moving in the trap, is a 
function of both position r{t) (due to the inhomogeneous 
field distribution) and time t (due to the rotation of the 
field vectors relative to the trap center), as discussed in 
Section [III The time dependent Hamiltonian that the 
molecule experiences is calculated by assuming a position 
in time given by r{t) = {vt — a;init)£ + by, which corre- 
sponds to the molecule moving with a constant velocity 
V and approaching the trap center with a minimum dis- 
tance b. The variable Xi^n must be chosen such that the 
initial position r(0) = —XiT^itX + hy is sufficiently far from 
the center such that at r(0) the adiabaticity condition is 
still satisfied. 

The initial state of the molecule could be chosen as a 
low-field-seeking eigenstate of the instantaneous Hamil- 
tonian at the initial position and time. If the electric 
field vector is rotating, however, the wavefunction of the 
molecule will immediately accumulate amplitude in other 
quantum states, even if the magnitude of the electric field 
is constant. The amplitude that appears in other states 
will only be negligible if the energy splitting between 
them and the initial state is much larger than the ro- 
tation frequency. Alternatively, one can choose an initial 
state that is a stationary state in the rotating system, as 
will be discussed below. The calculation of the transition 
probability as the molecule flies past the field minimum 
can then be started at lower fields, i.e. with a smaller 
value of Xinit- 

In the current system, the magnetic field vector is ori- 
ented along the z axis and the electric field vector lies in 
the xy plane. If the angle of the electric field vector with 
respect to the +x axis (toward the -\-y axis) is given by 
(/)', the Hamiltonian of the molecule can be computed by 
rotating the physical system around the z axis through 
an angle —(/>', operating on it with a Hamiltonian H' 
which corresponds to a system with the same magnetic 
field vector and in which the electric field has the same 
magnitude but is directed along the -\-x axis, and then 
rotating the physical system back through an angle 0'. 
In operator form, this can be written as 

H = e-'^-'t'' H'e'^-'t'' (5) 

where Fz is the total angular momentum of the molecule 
along the z axis. Using this form of H, equation ^ can 
be rewritten as 

where ipQ = e^^^'l' -0. This equation has the same form 
as equation ^ and can be solved in the same way as the 
time independent Schrodinger equation if the magnetic 
field is constant, the electric field strength is constant, 
and the electric field vector rotates at a constant fre- 
quency. The eigenvalues that result are called "quasiener- 
gies" and Hq = H' — F^H^ is the quasienergy Hamil- 



tonian [l^ . [TtI . A similar approach has been used in the 
recent work by Wall et al. Q. 

In section II it was shown that the angle of the electric 
field in the xy plane (j)' = (j)Q — (j) + it contains contribu- 
tions from both the angular coordinate of the molecule 
with respect to the trap center (p and a phase that re- 
sults from the constant rotation of the field vectors as 
the traps move over the chip 00 = 27r-^ t. Because 0o 
increases linearly in time, the contribution of 0o to the 
quasienergy Hamiltonian is time independent, having the 
form —h^Fz- This operator is diagonal in the basis 
sets used for both ^^CO and ^^CO, and in the case of 
^^CO, it is exactly equivalent to the Zeeman interaction 
for each energy eigenstate. As the rotation frequency 
^ is 3.75 MHz when the molecules are guided with 
300 m/s over the chip, its contribution to the quasienergy 
Hamiltonian is equivalent to a magnetic field of about 
—8 Gauss. It should be understood, however, that we are 
not dealing with a real magnetic field produced by the 
rotating electric field here; the effect of a rotating coordi- 
nate system merely produces a shift to the quasienergies 
that resembles the Zeeman Hamiltonian but that does 
not depend on the magnetic moment. In the case of 
^^CO, the rotation frequency is not equivalent to a Zee- 
man interaction since, as stated in the definition of /2, 
the gyromagnetic factors of the orbital angular momen- 
tum and the electron spin are different. 

For the subsequent calculations, the quasienergy eigen- 
vector instead of the normal Hamiltonian eigenvector is 
chosen as the initial state. After choosing an initial state 
i at the position r(0) = — Xiniti + by, the quasienergy 
vector is propagated in time using equation © until the 
molecule reaches the position r(tfined) = a^initi -I- by. The 
final state is then expressed in terms of quasienergy eigen- 
vectors at the final position and time, and the probabil- 
ity of the molecule ending up in a state j, Pij(v, B,b), 
is calculated. For all calculations, the population is as- 
sumed to be initially distributed equally over all low- 
field-seeking levels, two for ^^CO and four for ^'^CO. In 
^^CO, the calculation of Pij{v, B,b) only needs to be 
carried out for the initial state i corresponding to the 
Mj = +1 low-field-seeking level, since the Alf = 
low-field-seeking level is completely stable against non- 
adiabatic transitions (see the Appendix). For ^^CO, none 
of the four low-field-seeking levels is stable against non- 
adiabatic transitions at all magnetic fields, so Pij{v, B, b) 
must be calculated for each initial level i. 

In the end, we are interested in the probability T{B) 
of a molecule remaining in a low-field-seeking state for 
the duration of its time in the microtrap, as this is the 
quantity that is measured in the experiment. For ^^CO 
molecules in the ideal case, i.e. when the electric field 
is perfectly perpendicular to the magnetic field and the 
traps move perfectly smoothly over the chip, the sur- 
vival probability for a single molecule in a single state 
is given by the product of its survival probabilities af- 
ter each individual encounter with the trap center. The 
overall transmission probability T{B) is then calculated 
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by averaging this over both low-field-seeking states and 
over N molecules. 



TiB) = ^ XI ^ n PMf=o,Mf=oi^^,k,B, hr,.k) + 

n=l L fc=l 
1 ^ 



fc=i 



(7) 



The total number of passes K of each molecule and the 
speed Vn,k and closest approach distance bn,k of the kth 
pass of the nth molecule are determined using simulations 
of the classical trajectory of a molecule in the trap, as 
described elsewhere 01 • Since the Alf = state is com- 
pletely stable, PMf=o^Mf=oi'^-'n.k,B,bn,k) = 1, and the 
transmission probability T{B) in the ideal case can never 
be less than 1/2. For "CO, calculating T{B) IS some- 
what more complicated, since the F = 3/2, Mp = -1-3/2 
low-field-seeking level and the F = 3/2, = -1/2 low- 
field-seeking level are coupled, as are the F = 3/2, Mp = 
+1/2 and F = 3/2, = -3/2 levels. During each en- 
counter with the trap center, a molecule in a particular 
low-field-seeking state can transition not only to a non- 
trappable state but also to one other low-field-seeking 
state. To calculate the transmission probability of a 
single ^^CO molecule, the population in each low- field- 
seeking state after an encounter with the trap center is 
computed based on the population distribution before 
the encounter, and the total population still in a low- 
field-seeking state after the last pass is recorded. As in 
^^CO, T{B) is obtained by averaging the result of this 
calculation over a large number of molecules. 

To accurately describe the experimental data, the the- 
oretical calculations must be extended to include the jit- 
tering motion of the traps and the nonperpendicularity 
of the electric and magnetic fields. The jittering mo- 
tion is accounted for by including the full motion of the 
center of the trap as shown in Figure [3] in the calcu- 
lation of the transition probability Pi,j{v, B,b,t, (f>). As 
in the ideal case, the transition probability depends on 
speed u, magnetic field B, and closest approach distance 
b (although this distance is now defined relative to the 
average position of the minimum instead of the actual 
position). Additionally, the transition probability now 
also depends on the time t at which the molecule arrives 
relative to the jittering cycle and the direction </> from 
which it comes. If the electric and magnetic fields are 
not exactly perpendicular, low-field-seeking states that 
are normally decoupled can mix. In ^^CO, the Alf = 
and Mj = +1 states, which converge asymptotically at 
high electric fields, then become coupled. As a result, the 
population can partially redistribute between these two 
levels while the molecule is in a region of high electric 
field between successive encounters with the trap center. 
To account for this effect, it is assumed in the calcula- 
tions that, after each encounter with the trap center, a 
molecule's population n in each of these two low-field- 



seeking states is redistributed such that 

"Mf =+1 = (1 - rn)nj^.jB^^^ + muj^jB^f^ (8) 
'^Mf =0 = (1 - m)nMf=o + ™iMf=+i (9) 

where n' is the new population distribution. The pa- 
rameter m describes the degree of the redistribution; at 
the extremes, a value of indicates that no redistribu- 
tion occurs while a value of 1/2 corresponds to complete 
redistribution. Its exact value is difficult to predict and 
should actually depend on the trajectory of the molecule. 
For simplicity, m is determined by fitting it to the data; 
note that this is the only fitting parameter used. For 
^■^CO, remixing can occur at high electric fields between 
the F = 3/2, Mj? = -3/2 and the F = 3/2, Mj? = -1/2 
levels and also between the F = 3/2, Mp = +1/2 and 
the F = 3/2, = +3/2 levels. The remixing coeffi- 
cient m can be different for each of these pairs of levels, 
and thus for ^"^CO, two fitting parameters are necessary. 



IV. EXPERIMENTAL RESULTS 

To measure the effectiveness of the magnetic field 
at suppressing non-adiabatic losses, metastable CO 
molecules are guided over the full length of the chip at 
a constant velocity of 300 m/s and are subsequently de- 
tected using laser ionization. Measurements are carried 
out for both positive and negative magnetic fields, where 
the direction of positive magnetic field coincides with the 
-|-z-axis (see Figure [T]). To compensate for long term 
drifts in the intensity of the molecular beam, the parent 
ion signal is measured with the magnetic field on and with 
the magnetic field off, and the ratio between these two 
measurements is recorded. In Figure[Sl the thus recorded 
relative number of ^^CO (upper panel) and ^^CO (lower 
panel) molecules that are guided over the chip is shown 
as a function of the applied magnetic field. Both the 
measurements with the "standard waveforms" and with 
the "improved waveforms" are shown in the upper panel. 

It is clear from the data shown in Figure [S] that, as the 
magnetic field strength increases, the number of ^^CO 
molecules reaching the detector increases. This is as ex- 
pected because the splitting between the low-field-seeking 
levels of ^^CO and the non-trappable level increases with 
the magnetic field strength, thereby increasingly sup- 
pressing non-adiabatic losses. The measurements with 
the "improved waveforms" show that at magnetic fields 
more negative than —40 Gauss and more positive than 
-1-50 Gauss, the number of guided molecules becomes con- 
stant, indicating that all non-adiabatic losses are sup- 
pressed under these conditions. The data have been 
scaled vertically such that the saturation observed at high 
magnetic field strengths corresponds to a transmission of 
unity. The ratio between the signal at high magnetic 
field and low magnetic field is smaller than when the 
"standard waveforms" are used, implying that the "im- 
proved waveforms" also reduce the losses without a mag- 
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magnetic field [Gauss] 



FIG. 5: Measured signal of ^^CO (upper panel) and "CO 
(lower panel) molecules guided over the chip at 300 m/s as 
a function of the magnetic field. The vertical scaling of the 
individual data sets is described in the text. Measurements 
have been performed using the "standard waveforms" (filled 
circles) for both ^^CO and ^^CO; for ^^CO, measurements 
with the "improved waveforms" (open circles) are shown as 
well. The simulated curves of the transmission probability 
T{B) are shown as solid curves that overlay the data points. 
The thick solid curves in the upper and in the lower panel 
result from the theoretical model for the ideal case. The red 
dashed curves show the prediction of the theoretical model for 
the survival probability of molecules decelerated from 300 m/s 
to zero velocity in 250 ps. 



netic field. The number of CO molecules reaching the 
ionization detector, on the other hand, is seen to decrease 
as a function of magnetic field. The vertical scale in this 
case is based on the results of the theoretical calculations 
at low magnetic field strengths {vide infra). The lower 
panels of Figure 2] show that the Mp = +1/2 level of 
the non-trappable F = 1/2 state increases in energy in a 
magnetic field, while the = -3/2 and AI^ = -1/2 
levels of the low-field-seeking F = 3/2 state are at the 
same time lowered in energy. This reduction of the energy 
gap between the low-field-seeking and non-trappable lev- 
els enhances the non-adiabatic losses. Both the ^^CO and 
^■^CO data are slightly asymmetric for positive and neg- 
ative magnetic fields. The data for ^^CO indeed seem to 
be symmetric around a magnetic field of about +8 Gauss, 
as expected from the theoretical model. 

The results of the theoretical models are shown as solid 
curves in Figure [S] as well. The thick solid curves in the 



upper and in the lower panel result from the theoreti- 
cal model for the ideal case, i.e. when the jittering of 
the traps would be absent and no remixing would occur 
between decoupled states. The width of the theoreti- 
cally predicted transmission minimum for ^^CO around 
4-8 Gauss depends sensitively on the relative velocity be- 
tween the CO molecule and the center of the trap as 
they pass, and gets larger with increasing velocities. In 
the case of ^"^CO, two narrow transmission minima are 
expected around —65 Gauss and +90 Gauss, correspond- 
ing to fields at which the the F = 1/2, = +1/2 and 
F = 3/2, Mp = —3/2 levels cross, and between these two 
minima, no losses are expected. As in ^^CO, the width of 
the transmission minima increases as the relative velocity 
between the molecules and the trap center increases. It 
is clear from the comparison of these theoretical curves 
with the experimental data, that the observed measure- 
ments can not be quantitatively explained if the traps are 
assumed to move smoothly over the chip; the jittering of 
the traps must be taken into account. 

The transmission probabilities calculated for the case 
in which the jittering is explicitly taken into account 
are seen to almost quantitatively agree with the mea- 
surements. In particular the theoretical curves repro- 
duce the asymmetry between the intensity of the guided 
^^CO molecules at positive and negative magnetic fields 
as well as the narrowing of the profiles when the wave- 
forms are improved. In the calculations for ^^CO, a par- 
tial redistribution after each pass of 18% (m = 0.09) 
has been assumed for the "standard waveforms" and 12% 
(jTi = 0.06) for the "improved waveforms". It can not be 
excluded that the jittering of the traps was still slightly 
more severe during the actual experiments than shown 
for the "standard waveforms" in Figure [3l which would 
explain the experimentally observed additional broaden- 
ing for that case. In the ^■^CO calculations, the best 
agreement with the experimental data was found when 
assuming no remixing between the Mp = +3/2 and the 
Mj? = +1/2 levels, and 30% remixing (m = 0.15) be- 
tween the A/f = -1/2 and the M^ = -3/2 level. The 
maximum transmission at about +10 Gauss is not sensi- 
tive to the remixing coefficients, since the transition prob- 
ability for each of the low-field-seeking states is about the 
same in this region. It can be reliably inferred from the 
theoretical calculations that, even at low magnetic fields, 
about 1/3 of the ^■^CO molecules are lost to non-adiabatic 
transitions while being guided over the chip. Based on 
this, the ^■^CO data shown in the lower panel of Figure[5] 
have been scaled vertically such that the transmission at 
zero magnetic field is 2/3. 

The theoretical model used to explain the guiding data 
can also be applied to predict the non-adiabatic losses 
that are expected to occur during linear deceleration. 
The red dashed curves in the upper (lower) panel of 
Figure [5] show the survival probability of ^^CO (^'^CO) 
molecules decelerated from 300 m/s to zero velocity in 
250 lis. In these calculations, it is assumed that the jit- 
tering motion at 300 m/s is that of the "standard wave- 
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forms". The velocity of the jittering motion is assumed 
to be proportional to the frequency of the applied wave- 
forms while the latter is reduced from 2.5 MHz to zero. 
For ^^CO at low magnetic fields, the survival probability 
during deceleration is only 1/4 of the survival probabil- 
ity of ^^CO guided at a constant velocity of 300 m/s. 
The magnetic field needed to suppress losses is smaller, 
however. While the survival probability for guiding is 
symmetric around a magnetic field of +8 Gauss, the sym- 
metry point for deceleration is shifted closer to zero field, 
to +4 Gauss, due to the lower rotation frequency of the 
electric field vectors in the trap at lower velocities. For 
^■^CO, the model predicts that the transmission proba- 
bility during deceleration is larger than for guiding at all 
magnetic field strengths. 

The differences in transmission probability between 
guiding and deceleration can result from various effects 
that either enhance or suppress losses as the deceleration 
of the trap increases. The smaller spatial acceptance of a 
strongly accelerated trap results in a larger fraction of the 
trapped cloud being in the jittering region at any given 
time, which enhances the losses during deceleration 0|- 
Losses are also enhanced due to the molecules spending a 
longer time on the chip. On the other hand, since the av- 
erage velocity of a decelerating trap is lower than that of a 
trap at constant velocity, the velocity of the jittering mo- 
tion is reduced, suppressing non-adiabatic losses. While 
it is difficult to predict through simple arguments the 
relative importance of these effects, the outcome of the 
calculations is corroborated by previous measurements 
at zero magnetic field, in which it was shown that ^'^GO 
molecules can be decelerated to a standstill while ^^GO 
molecules are rapidly lost with increasing deceleration [l| . 



V. CONCLUSIONS 

In this paper we have studied the losses due to non- 
adiabatic transitions in metastable GO molecules — 
laser-prepared in the upper A-doublet component of the 
J = 1 level in the a^Hi, v = state — guided at a 
constant velocity in microtraps over a chip. Transitions 
between levels in which the molecules are trapped and 
levels in which the molecules are not trapped can be sup- 
pressed (enhanced) when the energetic splitting between 
these levels is increased (decreased) by the application of 
a static magnetic field. For a quantitative understand- 
ing of this effect, the energy level structure of ^^GO and 
^■^GO molecules in combined magnetic and electric fields 
has been analyzed in detail. When the GO molecules are 
guided over the chip, they are in an electric field that 
rotates with a constant frequency; the direction of the 
externally applied magnetic field is perpendicular to the 
plane of the electric field. The probability with which ei- 
ther ^^GO or ^"^CO molecules are transmitted over the 
chip, i.e. the probability that the molecules stay in a 
trapped level for the complete duration of the flight over 
the chip, has been measured as a function of the mag- 



netic field. The observed transmission probability can be 
quantitatively explained. 

To reduce trap losses in future experiments, it will 
be important to improve the applied waveforms. This 
will not only reduce losses due to non-adiabatic transi- 
tions caused by the jittering, but it will also reduce losses 
due to mechanical heating. Mechanical losses might also 
have been present in the current experiment, but because 
the measurements always compared the guiding efficiency 
with magnetic field on and off, we have not been sensi- 
tive to these losses. Alternatively, it might be possible 
to avoid the need for improved waveforms by moving the 
minimum on an orbit that is much larger than the am- 
plitude of the jittering motion, creating a large region 
around the effective trap center through which the mini- 
mum never passes. Such a trap, known as a time orbiting 
potential (TOP) trap, prevents non-adiabatic losses but 
is much shallower than a static trap 

The intrinsic difficulties with making the waveforms 
required for the experiments discussed here should be 
stressed; with present day technology these waveforms 
can hardly be made better than we have them now, in 
particular because, in order to bring molecules to a stand- 
still, we want to be able to rapidly chirp the frequency 
down from 2.5 MHz to zero. While the LG filter used 
to produce the "improved waveforms" reduces the total 
harmonic distortion of the amplitudes from 7% to 3%, 
it also makes producing a constant amplitude frequency 
chirp more complicated. We are nevertheless optimistic 
that the jittering can be reduced by another factor of 
two to three relative to the best waveforms that we have 
used so far. In the case of ^^GO, for instance, a mag- 
netic field of 10 Gauss, applied in the right direction, 
would then already completely avoid losses due to non- 
adiabatic transitions. With the present waveforms, trap 
losses can only be avoided when the applied magnetic 
fields are made sufficiently high. One should realize that 
there is an upper limit to these fields, however, as at some 
point transitions to the lower A-doublet components can 
be induced, opening up a new loss-channel. 

The extreme sensitivity to the details of the applied 
voltages results from the fact that the electric field min- 
ima above the chip originate from the vectorial cancel- 
lation of rather large electric field terms. Design studies 
are in progress to find an electrode geometry that is less 
sensitive to imperfections in the applied waveforms. A 
modified electrode geometry is also required to avoid trap 
losses at the ends of the tubular traps. Although the ends 
are closed in the present geometry by the fringe fields of 
adjacent electrodes, the electric field near the ends has 
components along the long axis of the trap, presumably 
leading to non-adiabatic losses even in the presence of 
the offset magnetic field. 
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Appendix: The a^IIi, v = 0, J = 1 Hamiltonian 

In this Appendix, the formalism that has been used to 
calculate the energies of the M -components of the J = 1 
level in the a^Hi, v = state of both ^^CO and ^^CO in 
combined, but mutually orthogonal, static electric and 
magnetic fields is presented. In the coordinate system 
used here, the magnetic field vector is oriented along the 
z axis and the electric field vector is in the xy plane. In 
this case, the molecular Hamiltonian is invariant under 
reflection in the xy-plane. It is thus possible to sepa- 
rate the basis states into two uncoupled sets, consisting 
of wavefunctions that are either symmetric or antisym- 
metric under reflection in the xy-plane, thereby reducing 
the computational complexity. The magnetic fleld can 
only couple states of the same parity and the same Mp 
quantum number, where Mp is the projection of the to- 
tal angular momentum including nuclear spin along the 
+z axis. In the case of ^^CO, there is no hyperfine inter- 
action and thus F = J and Mp = Mj . As the electric 
field vector lies in the plane perpendicular to the quan- 
tization axis it can only couple states of opposite parity 
with Mp differing by ±1. The two resulting sets of un- 
coupled basis states for ^^CO are given by: 

1. Mf = -l,± 

2. Mf = 0,T 

3. Mf = 1;± 

The + and — sign at the end describe the parity of the 
basis state. All states with the upper (lower) sign belong 
to one set. For ^^CO there are two sets containing six 
basis states each: 



1. 


F 


= 3/2, MP = 


-3/2, ± 


2. 


F 


= l/2,i\/J? = 


-1/2,T 


3. 


F 


= 3/2, A/# = 


-1/2,T 



4. F = 1/2, Aff = 1/2, ± 

5. F = 3/2, M^ = l/2,± 

6. F = 3/2, M^ = 3/2, T 

Again, all states with the upper (lower) parity belong to 
one set. 

Based on this formalism and using the zero-field spec- 
troscopicparameters and matrix elements given in ref- 
erences |8l4l5l|. the corresponding Hamiltonian matrices 
can be calculated. Without loss of generality, the elec- 
tric field vector is taken to be oriented along the x axis, 
i.e. E = Ex. The Hamiltonian matrices for other ori- 
entations of the electric field vector in the xy plane can 
be computed using the unitary transformation given in 
equation ([5]); the matrices given here correspond to H' 
in this equation. The origin of the energy scale for each 
isotopologue is defined to be the lowest energy field-free 
state in the upper A-doublet component, as shown in 
Figure U) 

For ^^CO, the matrices iJupper and Hiowa for the set of 
basis states with the upper and lower parity, respectively, 
are given by: 



Ri 


s 





S 


-A 


s 





S 


Ri 



and 

/-A-i?2 S* \ 

Slower =\ S S \ (A.2) 

\ S -A + R2J 

where A = 394.066 MHz 

Ri = {l,+\Hz\l,+) = 0.3332^lBB 
R2 = {l,-\Hz\l,-) =0M06fiBB 
S = (0, +\Hs\l, -) = 0.3513fiEE, 

and the bra and ket vectors have the form jAfJ^, parity). 
It is clear from these matrices, that the energy level la- 
beled as M f=0 in the case of ^^CO (upper right panel of 
Figure 2]) is only directly coupled to the Mj=ztl levels 
of the lower A-doublet component and that there is no 
direct coupling to the nearby Mf=±l levels of the up- 
per A-doublet component. Provided that the electric and 
magnetic fields are exactly perpendicular, this Mj=0 
level is therefore stable against non-adiabatic transitions. 

For ^■^CO, the corresponding matrices for the sets of 
basis states with the upper and lower parity are given 
by: 
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and 



ppcr 



1 — 


^/3S'3 


#54 








\ 




— i?3 — i?2 


-R4 


Si 


-53 





— D4 


i?4 


—£'2 — -Re 


S2 


*S'4 








Si 


52 


Ri 


R3 


-V3S2 





-S3 


54 


R3 


El + -R5 


2 











-V352 


2 '-'4 


-£2 + 3i?6/ 



(A.3) 



-f^lowcr — 



/- 


--E2 ^ S-Re 


V352 


#54 








\ 




V3S2 


~Ri 


R3 


Si 


-52 







#54 


R3 


El — -R5 


S3 


54 










Si 


53 


~E3 + R2 


-R4 


-V353 







-52 


54 


-R4 


—E2 + -Re 


■#54 













-V353 


vie-, 

2 "^4 


El + 3rJ 



(A.4) 



where Ei = 58.412 MHz 
E2 = 309.340 MHz 
E3 = 346.346 MHz 

Ri = (1/2, 1/2, +\Hz\l/2, 1/2, +) = 0.2264/xbS 
i?,2 = (1/2, 1/2, -|iJz|l/2, 1/2, -) = 0.2309Aii3B 
R3 = (1/2, 1/2, +\Hz\3/2, 1/2, +) = 0.1600Aii3B 
Ra = (1/2, 1/2, -|iJz|3/2, 1/2, -) = 0.1635AiBB 
i?5 = (3/2, 1/2, +\Hz\3/2, 1/2, +) = 0.1132/iBB 
i?e = (3/2, 1/2, -\Hz\3/2, 1/2, -) = O.llSS/isB 



51 = (1/2, -1/2, +|ffs|l/2, 1/2, -) = 0.3313AiB£; 

52 = (3/2, -1/2, -|if5|l/2, 1/2, +) = 0.n72^lEE 

53 = (3/2, -1/2, +\Hs\l/2, 1/2, -) = 0.1172AiB£; 

54 - (3/2,-1/2, +11^513/2, 1/2, -) = 0.3313^lEE, 
and the basis vectors have the form \F, Mj?, parity). 

If the magnetic field is not perpendicular to the electric 
field, additional non-zero matrix elements will appear in 
the Hamiltonian that couple the states of -ffuppor and 

-Slower • 
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